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The  work  reported  herein  was  performed  at  the  Metals  Behavior 
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under  Contract  F33615-77-C-5003  with  the  Southeastern  Center  for 
Electrical  Engineering  Education.  The  work  was  conducted  by  Dr.  Donald 
C.  Stouffer,  on  leave  from  the  University  of  Cincinnati  and  Or.  Sol  R. 
Bodner,  on  leave  from  the  Technion-Israel  Institute  of  Technology  in 
Haifa.  Air  Force  administrative  direction  and  technical  support  was 
provided  by  Dr.  Theodore  Nicholas,  AFML/LLN. 
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SECTION  I 
INTRODUCTION 

An  unresolved  problem  in  the  characterization  of  plastic  deformation 
of  metals  is  the  formulation  of  hardening  laws  that  would  realistically 
represent  hardening  properties  for  completely  general  loading  histories. 

Most  of  the  proposed  theories  follow  classical  plasticity  so  the  problem 
becomes  one  of  predicting  the  size  and  shape  of  yield  surfaces  subsequent 
to  loading  in  the  plastic  range.  These  laws  are  generally  some  combination 
of  the  isotropic  and  kinematic  hardening  models  and  reduce  to  either  of 
those  as  a special  case.  A feature  common  to  many  of  these  hardening 
laws  is  a "back  stress"  term  which  corresponds  to  a translation  of  the 
origin  in  the  kinematic  model.  Examples  of  these  hardening  formulations 
have  been  discussed  recently  in  the  proceedings  of  several  conferences 
(References  1,  2,  3),  and  others  have  since  been  proposed  (References  4,  5). 

For  constitutive  equations  that  do  not  require  a yield  criterion, 
the  corresponding  problem  is  the  determination  of  suitable  evolutionary 
equations  for  the  internal  state  variables.  These  variables  appear  in 
the  equations  for  plastic  deformation  rate  and  represent  the  load  history 
dependent  worked  state  of  the  material  with  respect  to  plastic  flow. 

Most  of  the  earlier  work  in  the  development  of  such  constitutive  equations 
employed  one  or  two  scalar  functions  of  deformation  history  for  the 
inelastic  state  variables.  As  an  example,  Bodner  and  Partom  (References 
6,  7),  used  a single  scalar  function  of  plastic  work  as  the  inelastic 
state  variable  for  isotropic  hardening  conditions,  e.g.  uniaxial  stress 
with  no  Bauschinger  effect.  In  this  theory  the  plastic  strain  rate  is 
assumed  to  exist  at  all  stress  levels  so  the  "back  stress”  concept  does 
not  appear  to  be  particularly  useful  since  it  would  correspond  to  a non- 
zero plastic  deformation  rate  at  zero  applied  stress. 

The  main  problem  in  generalizing  the  constitutive  equations  to 
multidimensional  stress  states  is  the  development  of  suitable  evolu- 
tionary equations  for  the  inelastic  state  variables.  These  equations 
are  required  to  account  for  the  loading  induced  changes  in  resistance 
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to  plastic  flow  which  vary  in  orientation  as  well  as  with  the  sign  of 
the  applied  stress  (or  plastic  deformation  rate).  The  material  therefore 
develops  both  orientational  anisotropy  and  a multidimensional  Baushinger 
effect  which  depends  upon  the  complete  loading  history.  Under  the  most 
general  stress  state  with  six  stress  components,  a minimum  of  12  inelas- 
tic state  variables  (components)  would  be  required  to  characterize  the 
material  state  with  respect  to  plastic  deformation;  this  assumes  a 
single  variable  would  apply  for  each  orientation  and  direction  of  the 
plastic  rate  of  deformation  component. 

The  purpose  of  this  report  is  to  develop  a proper  mathematical 
formulation  for  the  evolutionary  equations  required  to  represent  the 
inelastic  state  of  material  subjected  to  general  loading  histories  at 
large  deformations.  The  procedure  is  motivated  by  earlier  work  of 
Bodner  (Reference  8)  and  Onat  (Reference  9)  and  the  constitutive  equations 
of  Bodner  and  Partom  (References  6,  7,  10,  11)  but  the  development  is 
completely  general  and  could  be  applied  to  any  set  of  constitutive 
equations  that  employ  inelastic  state  variables.  A consequence  of  this 
development  is  that  a procedure  is  obtained  that  could  account  for 
induced  plastic  anisotropy  for  arbitrary  loading  histories.  The  limited 
comparisons  possible  with  experimental  data  indicate  that  the  proposed 
procedure  is  consistent  with  these  observations. 
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SECTION  II 

A CONSTITUTIVE  MODEL  FOR  ISOTROPIC  DEFORMATIONS 


Let  d..  represent  the  rate  of  deformation  in  an  Eulerian  coordinate 

* J 

system  at  any  time  t.  The  rate  of  deformation,  defined  as  the  symmetric 
part  of  the  velocity  gradient,  is  written  as 


d 


ij 


(1) 


the  composition  of  an  elastic  rate  of  deformation,  d®,,  and  an  inelastic 

n w 

(or  plastic)  rate  of  deformation  dP,  . The  elastic  component  is  fully 

• J 

recoverable  whereas  the  inelastic  component  is  nonreversible.  The 
decomposition  (Equation  1)  is  a very  fundamental  part  of  the  current 
development  since  we  assume  that  both  components  are  non-zero  under  the 
action  of  all  non-zero  deviatoric  stress  components.  A continuous 
representation  of  plastic  flow  is  thereby  obtained  without  the  use  of  a 
yield  surface  or  loading  and  unloading  conditions. 

This  study  is  directed  toward  developing  a representation  for 
anisotropic  plastic  flow  at  large  strains  where  the  plastic  anisotropy 
could  be  initial,  deformation  induced,  or  both.  Elastic  strains  for 
metals  are  generally  small  even  with  large  plastic  deformations  so  it 
would  be  generally  appropriate  to  use  an  isotropic  elastic  formulation 
if  the  material  were  initially  isotropic.  It  is  also  well  established 
that  large  plastic  deformations  have  small  effects  on  the  elastic  con- 
stants even  with  significant  induced  plastic  anisotropy.  The  mathe- 
matical treatment  of  the  elastic  deformation  rate  d*r . at  large  strains 
would  then  be  that  given  in  Reference  11.  Strong  elastic  anisotropy 
could  be  treated  by  classical  linear  anisotropic  elasticity  if  the 
elastic  strains  were  not  large.  The  case  of  strong  elastic  anisotropy 
and  large  elastic  strains  is  practically  rare  and  would  have  to  be 
treated  as  a special  problem. 
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The  governing  law  for  isotropic,  plastic  deformation  is  taken  in 
a form  similar  to  the  Prandtl -Reuss  flow  law 


ij 


4, 


X(J2,  zK,  T)SU 


(2) 


The  quantities  d^  and  S^j  are  the  deviatoric  components  of  the  rate  of 
plastic  deformation  and  stress  tensors.  Thus  Equation  2 implies  that 
isotropic  inelastic  deformation  is  isochorlc.  The  scalar  material 
function  \ is  assumed  to  be  a function  of  the  second  Invariant  of  the 
deviatoric  stress  tensor,  ^ the  absolute  temperature,  T,  and  a set  of 
internal  state  variables,  z^.  The  quantities  zK  are  used  to  describe 
the  state  of  the  material  microstructure  at  any  time  t;  thus,  they 
depend  on  the  history  of  the  deformation  up  to  the  current  time  t and 
are  given  by  an  evolution  equation  of  the  form 


•k*1'  ■ W V T) 


(3) 


To  begin  to  establish  a definite  structure  for  Equations  2 and  3, 
let  us  use  one  specific  component,  z,  of  the  vector  zK  as  a hardening 
term;  e.g.  an  increase  in  z corresponds  to  an  increase  in  the  resistance 
to  plastic  flow.  Thus  assume  the  hardening  rate,  z is  a monotonical ly 
increasing  function  of  the  plastic  work 


“p(t) 


■i  '* 


(T)dj  , (i)dT 


ij 


(4) 


Note  that  this  formulation  neglects  hardening  recovery  (softening). 
The  evolution  equation  can  be  written  as 


z(t)  = g(u>  ) 
P 


(5) 


which  can  be  expressed  in  the  variables  of  Equation  3 as  shown  in 
Section  V.  After  integrating,  Equation  5 gives 


z(t) 


s <- 


j i 
} 0 


z(x)dT 


(6) 


where  z is  the  initial  hardness  of  the  material. 
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Under  uniaxial  flow  conditions  the  material  develops  a directional 
characteristic  so  that  the  rate  of  deformation  is  different  for  loading 
in  tension  and  compression;  i.e.,  the  Bauschinger  effect.  This  effect 
can  be  expressed  by  introducing  two  hardness  variables,  z+(t)  and  z"(t) 
that  are  used  for  positive  (tensile)  and  negative  (compressive)  rates 
of  plastic  deformation,  respectively.  Mathematically,  then 

z(t)  = z+(t)  11  (djj)  + z"(t)  H (-d'n)  (7) 

where  d^  is  the  uniaxial  plastic  deformation  rate  component  and  H is 
the  Heaviside  unit  step  function.  Notice  that  using  Equation  7 in 
Equation  4 does  not  produce  a discontinuity  in  d^  since  the  jump 
from  z+  to  z',  or  vis-a-vis,  occurs  when  d^  is  zero. 

The  values  of  z+(t)  and  z’(t)  can  be  determined  from  Equation  6 once 
we  develop  a representation  for  z+  and  z , respectively.  Recall  z is 
defined  as  the  total  rate  of  hardening  due  to  an  axial  flow  d^ . 

Following  the  concepts  of  Reference  7 assume  d^  is  positive  and  let  qz 
and  (l-q)z  correspond  to  the  isotropic  and  directional  hardening 
fractions*  so  that  z+(t)  = q z + (l-q)z  = z(t)  and  z’(t)  = q z - (l-q)z  = 
(2q-l)z  . However  if  the  dP^  is  negative  rather  than  positive,  then 
the  hardening  in  tension  and  compression  is  given  by  (2q-l)z  and  z, 
respectively.  Thus  we  can  write 


• + 


z (t)  = q z + (l-q)z 


d?, (t) 


'11 


11 


(Ol 


and 


z (t)  = q z - (l-q)z 


dn(T) 

I di1(T) I 


(8) 


♦These  are  similar  to  isotropic  and  kinematic  models  of  yield  surface 
plasticity. 
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* ♦ 

If  g-1  then  t i and  the  model  corresponds  to  isotropic  hardening; 
however,  if  q-0  then  the  model  contains  only  directional  hardening. 

If  q-  0 then  the  model  predicts  softening.  Combining  Equations  6,  7, 
and  8 (lives 


z ( t ) - z + 


I 


q z ( t ) Ci  t + 


an(t) 

IdJ’jUf 


f ( l-q ) z ( t ) 

J 0 djjfl) 


d t 


(9) 


The  material  parameter  q could  be  approximated  as  a material  constant  but 
for  more  exact  representations  it  may  have  to  be  treated  as  a history 
dependent  internal  state  variable.  The  first  integral  represents  the 
total  isotropic  hardening  up  to  time  t,  whereas  the  total  directional 
hardening  used  to  characterize  the  Bauschinger  effect  is  given  by  the 
second  integral. 

In  the  following  Sections  the  above  representations  for  isotropic 
and  uniaxial  plastic  flow  will  be  extended  to  a full  three-dimensional 
anisotropic  theory. 


6 


4b 


AFML-TR-79-401 5 


SECTION  III 

CONSTITUTIVE  FRAMEWORK  FOR  ANISOTROPIC  PLASTIC  FLOW 


Consider  the  Prandtl-Reuss  theory  (Equation  2)  for  Isotropic  plastic 
flow.  These  laws  have  an  important  property;  namely  a given  deformation 
rate  component,  d^  for  example  is  parallel  to  the  corresponding  stress 
component  S^.  We  will  show  that  this  property  carries  over  to  all 
anisotropic  materials. 

Let  us  now  rewrite  the  above  theory  to  include  anisotropic  plastic 
flow.  The  most  general  anisotropic  flow  equation  can  be  written  as 


Xijkt  Sk£ 


(10) 


where  is  a fourth 

state  variables.  Using 


order  tensor  valued  function  of  Jg,  T,  and 
the  symmetry  of  d!?j  and  S^.  it  follows  that 


Xijk!l  X j IkJl  Yjtk 


(11) 


Mathematically  and  conceptually  it  is  convenient  to  replace  the 
deviatoric  stress  and  rate  of  deformation  tensors  in  Equation  10  by 

A 

six-dimensional  vectors  on  the  vector  space  i , ct=l,2,...6.  Thus,  let  us 
define 


D = D i and  T = T i 
_ a a _ a a 


(12) 


and  normalize  the  six-dimensional  vector  space  such  that 


3 3 2 

l l (dj  ) and 
1=1  j=l  3 


3 

l * 

i-1  J-l 


(13) 


I 


h 

N 


' 
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One  such  mapping,  satisfying  Equation  13  is  * 


Using  this  notation  the  plastic  work.  Equation  4 can  be  rewritten  as 


u> 

p 


f D(t) 
J O' 


T(t)  ill 


(15) 


The  components  of  each  of  the  two  vectors  T and  D are  linearly  independent, 
as  such  the  vectors  T and  D are  orthogonal  Cartesian  vectors.  If  we  let 
a represent  a three-dimensional  proper  orthogonal  transformation  such 
that  S'  - a S n^,  for  example,  then 


T*  * C T and  O'  -CD 


(16) 


where  C is  a six-dimensional  proper  orthogonal  transformation  that  can 
be  easily  constructed  from  the  components  of  a. 

Let  us  now  rewrite  Equation  10  on  the  six-dimension  of  vector  space  as 

? - i T or  Da  ■ W»  (17) 

*A  mapping  similar  to  this,  except  on  a five-dimensional  space,  is  given 
by  Ilyushin  and  Lensky  (Reference  12). 


i 

i 

i 
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where  A is  a six-dimensional  square  tensor.  Let  us  assume  that  T can  be 
determined  from  a potential  function.  Then  using  commutativity  of  the 
partial  derivative,  it  follows  that 


A 


:t 


afl 


(18) 


Since  A is  symmetric  and  the  components  are  real,  there  exist  real 


eigenvalues,  A 


mo, 


Thus  using  A 3 0 A Q , where  Q Is  a proper  orthogonal 


transformation.  Equation  17  becomes 

I)  - A T 


(19) 


where  A is  a 6 x 6 diagonal  tensor.  The  rate  of  deformation,  D,  and 


stress,  T,  are  vectors  that  use  the  eigenvectors  e of  A as  a basis. 


Equation  19  has  the  advantage  that  all  the  anisotropy  properties  can 


be  described  by  the  six  eigenvalues  A and  the  transformation  L)  so  that 


the  structure  of  the  equation  is  particularly  simple.  The  details  for 
determining  Q are  delayed  until  the  end  of  the  next  section. 


Next  we  assume  that  each  non-zero  component  of  A,  can  be  expressed 


as  a single  valued  scalar  function  of  the  scalars  J^,  T,  and  Z^,.  where 


Z<(,  * 0 for  ot  / b,  so  that  each  non-zero  component  Z^,  corresponds 


directly  to  a non-zero  component  of  Mathematically,  this  Is  written 

as 


A - f (.)  ,T,Z  ), 

Old  '2  drt 


no  stun  on 


(20) 


Since  \ are  diaqonal  elements  of  a second  order  tensor  in  the  e basis, 
oux  • a 

then  Equation  20  implies  that  Z must  also  be  diagonal  components 

of  a second  order  tensor  on  the  basis  e , thus 

ex 


C 7 C 


„T 


(21) 


In  the  isotropic  theory  relating  to  Equation  2 we  imposed  the 
condition  that  the  plastic  flow  is  Isochoric.  Kor  anisotropic  flows  this 


Af  ML  - TR- 79-4015 


condition  leads  to  considerable  mathematical  difficulties  as  well  as  a 
physically  inconsistent  theory.  Using  Equation  19  consider  the  sum 


+ D1  ‘ 


Si  s 


If  the  material  is  isotropic,  then  3 X ^ 3 X^  * X and  Equation  22 
becomes 


D1  * D2  * "3 


,<sn  * s22  * V 


for  every  choice  of  X.  Thus  the  anisotropic  flow  model  will  reduce  to 
classical  Incompressible  isotropic  plasticity.  However  the  volumetric 
change  + D,,  + will  not  be  Identically  zero  for  anisotropic  flows. 
This  appears  to  be  consistent  with  experimental  observations,  since  the 
plastic  Poisson  ratio  measured  during  plastic  deformation  Is  generally 
different  from  one  half  (Reference  15). 
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SECTION  IV 

STRUCTURE  OF  THE  STATE  VARIABLES 

To  establish  a representation  of  Z for  anisotropic  plastic  flow 
we  follow  the  concepts  underlying  Equation  9.  This  was  used  to  intro- 
duce a nonsymmetric  property  Into  the  plastic  flow  equation  so  that  the 
response  In  tension  and  compression,  or  shear  to  the  left  and  right,  are 

different.  To  maintain  this  characteristic  let  us  define  Z in  the 

aa 

following  manner 

z (t)  - z+  (t)H[n  (t)l  + z"  (t)Hl-D  (t ) ] (24) 

aa  au  a aa  a ' ' 

where  Z+  (t)  and  Z"  (t)  are  the  hardness  variables  for  positive  and 

aa  mot  r 

negative  plastic  deformation  rates. 

It  is  informative  to  consider  three  different  classes  of  mechanical 
properties . 

a)  A material  is  defined  as  orientational ly  isotropic  if 


(25a) 

(25b) 


for  all  time  and  deformation  histories.  The  quantity  I is  a six- 
dimensional identity  tensor, 


!1  a - 0 

0 a i 0 


i 

I 

Thus,  according  to  Equations  25  and  23  for  equal  values  of  stress,  T , 
the  deformation  rates  are  the  same  in  all  six  coordinate  directions  of  the 
same  sign.  Further  using  either  Equation  25a  or  25b  in  Equation  21  gives 

(Z+)  - C(z+I)cr  - Z+I  (27) 
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which  shows  that  Z (and  l)  are  Identical  In  all  coordinate  directions. 
Note  that  this  result  follows  since  Z transforms  as  a second  order 
diagonal  tensor. 

b)  A material  is  skew  symmetric  if 

C(t)  ■ za6(t)  (28) 

for  all  time  and  deformation  histories.  This  corresponds  to  no 
Bauschinger  effects  in  a uniaxial  deformation  history.  Using  Equations 
19,  20,  and  28  we  can  write,  for  example, 

f(J2.  T,  Z~  ) (-T  > - -f(J,,  T,  Z+  )T 

- an  aa  2 aa  aa 

which  is  skew  symmetry  of  the  material  function  f. 

c)  A material  is  defined  as  Isotropic  If  It  Is  both  orientational ly 
isotropic  and  skew-symmetric . 

Consider  now  the  internal  state  tensor  Z.  Since  the  only  non-zero 

A 

components  of  L are  on  the  diagonal,  Z.  are  eigenvalues  of  Z and  e are 

OKI  -v  01 

the  eigenvectors.  Thus  Z can  be  written  as  a six-dimensional  vector 
in  the  form 


Z(t)  - Z (t)  e (2 

- out  a 

Referring  to  the  structure  of  Equation  9,  let  us  rewrite  the  vector  Z 
in  the  same  general  form. 


Z(t)  - 9 + [ Zl(i)di  + f ?.A( i )di 

J i)  ~ J it  - 


o 

The  vector  Z is  used  to  specify  the  initial  hardness  in  each  of  the  six 

, I , a 

coordinate  directions.  The  vectors  Z and  Z correspond  to  the  rate 
of  fully  isotropic  and  directional  hardening,  respectively.  Note  that 
isotropic  hardening  is  characterized  by  the  first  integral  in  Equation  9 
whereas  the  second  integral  in  Equation  9 characterizes  orientat ional ly 
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anisotropic  and  non-skew  symmetric  hardening.  Thus  we  can  use  the  first 
integral  from  Equation  9 directly  into  Equation  30  to  obtain 

7.(t  1 « 7 + e | q z(Odi  + f 7 (O  dx  (31) 

a J 0 /O' 

where  2 ( x ) is  calculated  from  Equation  5 using  the  total  work  as  given 
by  Equation  lb.  Recall  q is  a material  parameter  describing  the  relative 
amount  of  fully  isotropic  hardening. 


Using  Equation  19  it  can  be  seen  that  there  is  a one  to  one  relation- 
ship between  a deviatoric  stress  component  and  the  corresponding  rate 
of  deformation  D^.  However,  it  is  not  expected  that  the  vectors  T and 
D are  parallel  since  the  components  of  Z,  and  hence  A,  are  not  neces- 
sarily equal.  Let  us  assume  that  the  rate  of  change  of  resistance  to 
plastic  flow  (hardness)  is  due  essentially  to  the  plastic  straining. 

That  is,  the  change  in  resistance  to  plastic  flow  is  in  the  direction 

• A 

of  the  rate  of  plastic  deformation  D;  so  that  Z and  D are  parallel. 
Expressed  mathematically 


2* 


(32) 


where  u are  the  "direction  cosines"  of  the  rate  of  deformation  vector. 
a 

Hence  we  can  write 


u 

it 


D 


(33) 


1 /2 

and  | D | = (D* D)  ' . This  formulation  gives  the  hardness  as  varying 

according  a cosine  rule  as  proposed  by  Bodner  (Reference  8).  Using 

A # 

again  the  second  integral  of  Equation  9 and  noting  that  |Z  | = (l-q)z(i) 
we  can  write 


7(t) 


. V‘>  r‘ 


q z(i) 


d i 


(l-q)z(l) 


+ 


U ( 1 ) d t 


(34) 


13 


I 
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Thus  wo  have  a specific  form  for  the  evolution  of  the  tensor  Z.  This 
representation  Is  based  on  the  main  assumption  that  the  rate  of 
anisotropic  hardening  vector  is  parallel  to  the  rate  of  plastic  defor- 
mation vector.  The  total  anisotropy  at  any  time  depends  on  the  integral 
of  Z(i)  or  the  history  ot  the  deformation.  Thus,  the  representation 
has  the  necessary  path  dependence. 

Returning  to  Equation  19,  it  is  necessary  to  know  the  transformation 
Q in  order  to  obtain  a diagonal  form  for  the  material  tensor  \.  From 
Equation  34  it  can  be  seen  that  Z is  written  relative  to  the  basis  e^, 

the  eigenvectors,  which  are  fixed  during  the  deformation  history  (0,t). 

o 

Thus  the  choice  of  Q depends  only  upon  the  initial  anisotropy  Z,  which 
is  assumed  to  be  a known  initial  condition.  Further  if  the  material 
is  initially  isotropic,  then  every  set  of  vectors  e are  eigenvectors 

CL 

and  the  vectors  D and  T can  be  chosen  such  that  the  components  are 
physically  convenient  for  the  shape  of  the  body. 

There  is  a limited  amount  of  experimental  evidence  available  for 
the  anisotropic  deformation  of  beryllium.  In  general,  Nicholas 
(Reference  13)  has  shown  that  a comprehensive  anisotropic  hardening 
theory  is  necessary  to  predict  the  experimental  results  for  beryllium. 

The  results  published  by  Lindholm,  Yeakly,  and  Davidson  (Reference  14) 
for  beryllium  in  biaxial  tension  showed  that  the  above  theory  is 
qualitatively  consistent  with  the  data  in  that  the  stress-strain  curves 
in  the  plastic  range  are  ordered  with  loading  history  according  to  the 
cosine  law  given  by  Equation  32. 


It 
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SECTION  V 

A SPECIFIC  REPRESENTATION 

The  purpose  of  this  Section  is  to  briefly  summarize  the  governing 
equations  for  deformation  Induced  anisotropic  plastic  flow  in  a form 
consistent  with  the  Bodner-Partom  model  for  isotropic  plastic  flows. 
Squaring  Equation  2 gives 

D2  ’ ‘2<J2-  V T)J2 


where  and  are  the  second  invariants  of  d?j  and  S^j,  respectively. 
Substituting  Equation  35  into  Equation  2 gives 


dP(J2,  zk,  T) 


This  result  is  reminiscent  of  the  Prandtl-Reuss  theory  except  that  no 
yield  surface  is  used  with  Equation  36;  i.e.,  is  continuous.  Using 
one  state  variable  and  motivated  by  the  equations  of  dislocation 
dynamics,  Bodner  and  Partom  selected  such  that 


dU  ■ °o  [-*  ^ 


where  the  constant  DQ  is  the  limiting  strain  rate  of  the  material  and  n, 
generally  a function  of  temperature,  relates  to  the  strain  rate  sensi- 
tivity. Since  z is  a single  internal  variable  that  controls  the  rate 
process  it  clearly  must  be  a measure  of  the  resistance  to  plastic  flow 
of  the  material.  Specifically  Equation  5,  is  chosen  in  the  form 


m(t1  - zn)  exp  [-m  u>^(t)]w^  = mlz,  - zlS^  d?j  (38) 


where  m,  z^,  and  Zq  are  constants  controlling  the  rate  and  level  of  the 
hardening  characteristics.  Equation  38  is  consistent  with  the  general 
form  for  z.  Equation  3. 
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The  representation  for  anisotropic  flows  can  be  completely  determined 
by  combining  the  above  isotropic  representation  with  the  results  in 
Sections  111  and  IV.  The  plastic  deformation  rate,  Equation  19  becomes 


where  D and  T are  the  six-dimensional  rate  of  deformation  and  deviatoric 

Ol  01 

stress  vectors  defined  by  the  mapping  (Equation  14).  The  components  of 
the  hardening  tensor  are  obtained  by  substituting  Equation  38  Into 
Equation  34. 

The  anisotropic  plastic  representation.  Equations  39,  38,  and  34, 
requires  the  use  of  five  material  constants  (q,n,m,z1 ,Zq)  that  must  be 
determined  experimentally  (Dq  is  a scale  factor).  It  is  useful  to  note 
that  these  material  parameters  can  all  be  determined  from  a one-dimensional 
experimental  program.  There  were  no  material  functions  introduced  by 
the  anisotropy  development.  This  leads  to  a specific  model  that  can  be 
readily  evaluated  for  different  materials. 

In  closing  it  is  important  to  note  that  the  constitutive  equations 
are  a nonlinear  coupled  system  that  can,  most  likely  be  solved  numerically. 
Numerical  finite  difference  and  finite  element  procedures  have  been 
developed  to  solve  boundary  value  problems  under  isotropic  hardening 
conditions.  It  is  expected  that  these  codes  can  be  modified  to  include 
the  anisotropic  model;  however,  to  date  no  numerical  work  has  been  done 
using  the  above  system  of  equations. 
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